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In this paper we discuss a Schwinger-Dyson [SD] approach for determining the time evolution of 
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. the unequal time correlation functions of a non-equilibrium classical field theory, where the classical 
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l/^ I system is described by an initial density matrix at time t = 0. We focus on A</)^ field theory in 1+1 

space time dimensions where we can perform exact numerical simulations by sampling an ensemble 

^ : 

■ of initial conditions specified by the initial density matrix. We discuss two approaches. The first, 

the bare vertex approximation [BVA] , is based on ignoring vertex corrections to the SD equations 
Q^l in the auxiliary field formalism relevant for 1/A^ expansions. The second approximation is a related 



approximation made to the SD equations of the original formulation in terms of (p alone. We com- 
pare these SD approximations as well as a Hartree approximation with exact numerical simulations. 



^ I We find that both approximations based on the SD equations yield good agreement with exact 

' numerical simulations and cure the late time oscillation problem of the Hartree approximation. We 

also discuss the relationship between the quantum and classical SD equations. 
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I. INTRODUCTION 



These past few years, there has been a concerted effort to find approximation schemes for 
time evolution problems that go beyond mean field theory, such as the Hartree|jl|, 0| or the 
leading order in approximations 0. Although these approximations are able to give a 
reasonable picture of the phase diagram of quantum field theories (see for example [Q), there 
are no true hard scatterings in mean field theory, so there is no mechanism for a system which 
is out of equilibrium to be driven back into equilibrium. Thus various phenomena found in 
these early simulations, such as production of disoriented chiral condensates [DCC's], and 
distortion of the low momentum spectra away from that of thermal equilibrium, might very 
well be modified when hard scatterings are taken into account. This is especially true in the 
phenomenological 0(4) sigma model for the chiral phase transition where the effective low 
energy couplings are large. In Relativistic Heavy Ion Collider [RHIC] physics applications 
it is important to know, given an expanding plasma, how the time scale for equilibration 
compares with the time scale for the expansion of the plasma in order to decide whether 
effects found in the mean field approximation persist following the chiral phase transition. 

Initial attempts to go beyond mean field theory, such as truncating the coupled hierarchy 
of n-point functions at the 4- or 6-point function level 0, as well as naive use of the 
expansion[0, |^, have suffered from serious drawbacks such as negative definite probability 
and/or secular behavior. These problems, although present in naive perturbation theory, 
were not present in the mean field approximation because the mean field approximation can 
be shown to be equivalent to a Hamiltonian dynamical system. Because of these failures, 
recent attention has been shifted to resummation methods based on Schwinger-Dyson [SD] 
equations resulting from keeping leading terms in the series for the two particle irreducible [2- 
PI] generating functional that occurs in the $ derivable approach of BaymQ, or equivalently 
the formalism of Cornwall, Jackiw, and Tomboulis [CJT] [T^, |ll|, |12|, |13|. In this paper we 
will restrict ourselves to looking at two proposed resummation methods based on the 2-PI 
approach. The first approach, is a direct use of the loop expansion for the 2-PI generating 
functional espoused by CJT. This approach has been used to study thermalization in the 
quantum version of the problem discussed here by Berges and Cox[|l^]. The second method 



is to first obtain the exact SD equations in the auxiliary field formalism |T5[, and then 
to make the assumption that one can ignore vertex corrections. This approximation also 
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can be obtained from the 2-PI formalism when we include propagators for the auxiliary 
field. This latter method we recently applied successfully to the quantum roll problem 
in 0+1 dimensions ||16|| , where we found that the approximation is a resummation of the 
next-to-leading order expansion, which cures the defects of the large- expansion in 
next-to-leading order. 

In this paper we will apply these approximations to a field theory problem where exact 
calculations can be carried out. Namely we will consider the time evolution of an initially 
Gaussian ensemble of classical fields in 1+1 dimensions. We will compare the results of the 
SD equation approximations with the exact simulation as well as with the Hartree approxi- 
mation. To make contact with earlier work in this area by Aarts, Bonini, and Wetterich^, 
which went beyond Hartree approximation by truncating the hierarchy of Green functions 
at the four point function level, we will calculate the same variables used by those authors 
as well as use their initial conditions. Since the quantum evolution in the high tempera- 
ture limit reduces to this classical field theory problem, by showing that the bare vertex 
approximation [BVA] gives reasonable quantitative agreement for the classical evolution we 
guarantee that the BVA gives reasonable results for the quantum evolution in the high 
temperature regime. 

The BVA approximation we espouse here is certainly not new. It was first introduced by 
Kraichnan||17|| as the "direct interaction approximation," which he used to study classical and 
quantum dynamical systems. What Kraichnan showed was that this approximation, like the 
Hartree, was realizable as a dynamical system and thus was free from the secularity problem 
as well as having a positive definite density matrix. A method of deriving these classical 
SD equations was obtained in the classic paper of Martin, Siggia and Rose [MSR] [ll8| . 
Although quantum field theorists are now familiar with the Schwinger-Keldysh closed time 



path [CTP] formalism for initial value problems [p!9|] , the related formalism of MSR is less 
familiar, and so we devote an appendix to deriving the exact SD equations for both classical 
as well as quantum dynamical systems and discuss the connection between the CTP and 
MSR formalisms. A more extensive dicussion of this connection has just been published!^ . 



This paper is organized as follows. In section II we discuss the classical field theory 
simulation we are using as a benchmark for various approximations. In section JTT| we review 
the Hartree approximation. In section |rV|, we discuss the BVA from various perspectives and 
compare the quantum and classical versions. In section |V|, we discuss the 2-PI approximation 
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used by Berges and Cox, and show that it is a re-expansion of the BVA. In section ^ we 
compare all three approximation schemes with the exact numerical solutions. We conclude 
in section |V11| . 

In addition, in the appendices we discuss both the CTP and MSR formalisms for a 
problem with two interacting fields and discuss the relationship between the quantum and 
classical SD equations. 

II. CLASSICAL SCALAR FIELD THEORY 
A. The Lagrangian 

We are interested in studying the classical scalar field theory in 1+1 dimensions described 
by the Lagrangian density: 

c = \[{^,<pf-{^.A?-^^'<P']-\<P\ (1) 

with the equation of motion: 

[d't-dl + i^'] 0(x,t) + ^0=^(x,t) = O. (2) 



In momentum space, we define: 



^'^ ' - Jkx 



<P{x,t)= I —<p,{t)e'>^\ (3) 



oo 



in which case (|^) becomes: 



^ J -oo 

where ul = k"^ + fi'^ . Because we are particularly interested in SD equations which are a re- 
summation of the large- iV expansion, we will also be using a second form of this Lagrangian: 

^ = ^[Wf-(5.0f-X0^]+f (f-/^^) • (5) 
This second Lagrangian leads to the equations of motion: 

[d^-dl + x{x,t)] 0(x,t) = O, (6) 

and the constraint ("gap") equation for x(x,t): 

X{x,t)=fi^ + ^(l)\x,t). (7) 



The Hamiltonian density is given by: 



n = - [{Mf + {d^f + ^ '^'l - f (f - /^') • (8) 



B. Initial ensemble 

Solutions of the equations of motion are found for the initial conditions 0) = 0(x) 
and 7r(x, 0) = Tt{x). If the system is in thermal equilibrium, the initial values (j){x) and tt{x) 
are taken from a canonical ensemble governed by a classical density distribution p[(/),7r], 
given by: 

p[0,7r] = Z-i(/?)exp{-/?i7[0,7r]}, 
Z{P) = 11 [[ ^^(^) exp{-/5i/[<^, vr]} . (9) 

X 

with /3 = 1/T. The ensemble average of a quantity y4[</), tt] is then defined by: 

(A[0,7r])=Tr{pA[0,7r]} 

^Z-\P)l[ /"/"d0(x)d7r(x)A[0,7r] exp{-/3i7[0,7r]}, (10) 

In order to make direct comparisons with the work of Aarts, Bonini and Wetterich[0, 
we choose an ensemble constructed from unperturbed Hamiltonian Hq rather than H in 
this paper. This is not a canonical equilibrium distribution for either the exact solution 
or the Hartree approximation, so it gives a non-trivial initial value problem for all our 
approximations. In classical 1+1 dimensional field theory, mass renormalization is finite 
so that one can use any arbitrary mass in our initial Hq] however, in order to compare 
our results with Aarts, et. ai, we take this mass to be the bare (unperturbed) mass fi. In 
dimensions higher than 1+1 one would need to use a renormalized mass parameter in Hq. 
So we choose for our initial density matrix: 

Po[0o, TTo] = ZQ^{f3o) exp{-/3oi^o[0o, TTq]} , (11) 

where Pq = I/Tq is to be regarded as a parameter. The unperturbed Hamiltonian is 

1 /'+°° 

i^o[0o, ^o] = 2 / [nlix) + {8^^))^ + /x^ <pl{x)] , (12) 
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with the equations of motion, 

[dl-dl + i?] 0o(x,t) = O, (13) 
where 7ro(x, t) = <po{x,t). Again, defining Fourier transforms by: 

— Mt)e''\ <x,t)= ^^,(t)e^'=-, 
solutions for the free particle case are given by: 

0ofc(i) = 0fc cos{uJkt) + (7rfc/a;fc) sin(u;fet) , 

T^Qkit) = TTk cos{uJkt) - (pk^k sin(u;fct) 



1 

7VT 



afc e-*^'^* - alfc e"^^' , (14) 



where cufc = \Jk? + /i^- Note that 0) = (^ol^^jO) and 7r(x, 0) = 7ro(s,0). So then the 
density matrix becomes: 

1 f f^"^ dk ^ 

Po [xk, Vk] = -^^^ I y 2^ + J ■ (15) 

where we have put ak = Xk + iyk- So for the Monte Carlo calculation, we select value Xk 
and i/k from the Gaussian distribution in Eq. (|l5|), and use these values to construct starting 



values for 00(2^5 0) and 7ro(x, 0). These functions are not smooth functions of x. From (|15|), 
we find that: 

(«X')o = (%'«fc)o = 27r(5(fc-fcVfc(/3o) , 

and where nk{Po) = l/{PoUJk)- Note that nk{Po) is the high temperature limit of the clas- 
sical Bose-Einstein occupation number distribution. Then, one show that ( 0(x, 0) )^ = 
( 7r(x, 0) )^ = 0, and that 



(o)4rd.(^=(.,o))„^rf!^ (17) 

L Jo J_A 27r cufc 2/i/?o 



<m = j dx(7r2(x,0))^ = j ^ —nk{(3,)ujk- ^ ^ ' 



Po 



[18) 
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where we have introduced cutoffs L and A in coordinate and momentum space integrals. 
We have then: 

2yU J_A 27r 7r2/i 
C. Correlation and Green functions 

We define correlation functions and Green functions by ensemble averages: 

F(x,x') = (0(x)0(x'))q, 
G'r(x, x') = — (t(x, x ) 0(t — t') , 

Ga{x,x') = +a{x,x')Q{t' -t) , (20) 

where (t(x, x') is the spectral function, given by the ensemble average of the classical Poisson 
bracket of with 0(x'): 

a(x,a;') = ({0(x),0(a;')}V (21) 

[In this section, we use x = (x, t) when needed.] Since po[0o,7ro] is Gaussian, the only 
nontrivial correlation and Green functions for the free particle case at t = are the two- 
point functions. These functions are discussed in detail in Ref. |^. Writing (I)q{x) in a 
Fourier expansion and using Eqs. (|16|), the unperturbed correlation and structure functions 
are given by: 

Fo(x, x') ^ { Mx) M^') )o = r ? ^^(^o) '"'^""'^'"'^^^ ' 

y_oo 27r ujk 

ao(x, X') ^ { {Mx), M^')} )o = - -^-M^-n] ^..(.-.0 _ (22) 

^ J-oo 27r iOk 

There are similar correlation and Green functions for the xi^) field. Even though this 
field is a constraint field, we treat it as if it is dynamic and then take a limit so that the free 
particle Green functions, which in this case we call Dor(x,x') and Doa{x,x'), become just 
delta functions. 
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D. Classical solutions on the lattice 



We introduce a lattice in coordinate and momentum space by putting and 



k ^ km, where: 



km = 27im/L 



n = 0, 1, 



m 



-N/2,-N/2 + l,... , N/2-1 



(23) 



where that L = Na and the momentum cutoff A = vr/a. We write (/>„(t) = (f){xn,t) and 
(pkmit) = (pmit)- Fourier transform relations become: 

N/2-1 N-l 



m=-N/2 n=0 

There are two strategies we can use to step out solutions. In coordinate space, we can 
introduce a time grid, ts = soq, s = 0, 1, . . . , and replace the differential operators for space 
and time by second order difference formulas: 



dU{x,t) 



^m,s+l 2 (f)m,s ~l~ 0m,s— 1 



So that Eq. (||) becomes: 



,2 



+ 



+ 



ao 



a 



{ 



>n+l,s 



+ 



/a^ 

/ 2 
^ ^3 



(24) 



}• 



(25) 



This method is called a "staggered leapfrog" approximation in the literature and is 
believed to be stable. The advantage of this coordinate space method on the lattice is that 
the interaction term, 0^(x,t), is easy to calculate. We construct an average of all Monte 
Carlo runs M for all values of x„ and compute: 

M N-l , „i 



j=l n=0 



(26) 



The use of the second order difference formula means that the dispersion relation is now 
given by: 



2 -^2 



km = — sin I 1 



/7rm\ 



(27) 
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So at t = 0, ( 0) )^ is still given by Eq. (|T^ but /(A) is now given by: 

We call this method "coordinate-space Monte Carlo." 

The second strategy, is to solve Eq. in momentum space, and use a fast Fourier 
transform to compute the convolution integral. Any time integrator can be used for the 
time step, but we employ the same second order finite difference formula as above. Here, 
there is no use of the second order difference formula, so the dispersion relation is given by: 

= km + , km= -J^ ■ (29) 



0) )^ is again given by Eq. (|T^) but /(A) is now given by: 

Af/2-l ^ , N/2-1 

^ '~ L ^ u;l~ [ 2A J ^ (nmf + (nfiN/2Ay 

m=~N/2 "I \ / ^^_^/2 ^ > ' \ I I 

- \ • (^°) 

for sufficiently large iV. 

We call this method "momentum-space Monte Carlo." Since our BVA codes use a global- 
in-time Chebyshev expansion method in momentum space, with a dispersion relation given 
by Eq. (|29|), we can compare directly the results of the BVA code to the classical momentum- 
space Monte Carlo calculations. 



III. HARTREE APPROXIMATION 



The Hartree approximation is a simple truncation scheme which is obtained by setting all 
correlation functions beyond the second one to zero. Since it has a simple interpretation in 
terms of a time-dependent variational approximation where the density matrix is assumed 
to be Gaussian, it automatically corresponds to a Hamiltonian dynamical system so it is 
free from unphysical behavior. Here we discuss the classical version of this approximation. 
We are interested in ensemble averages of the classical evolution. The equation for the 
expectation value 0c(a;) = )g is: 

{dl-dl^^^\ 0c(x) + ^(0^(a:))^ = O. (31) 
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[Here x = {x,t).] Since the Hartree approximation is equivalent to an initial Gaussian 
ensemble staying Gaussian, the expectation values obey the factorization condition 

where 

x') = { 0(x)0(x') )^ - ( 0(x) I { <t>ix') )^ . 

This equation must be supplemented by the equation for the expectation value of the fluc- 
tuation, which under the Gaussian hypothesis obeys: 

I - 5^ + /i^ + y { <j>l{x) + a;)} I G{x, x') = . (32) 

We will determine the 2-point function by assuming that we can write a mode decomposition 
for (f){x, t) of the form: 

— [a,/,(x)+a*/;(x)] (33) 

where the complete orthonormal set of mode functions fk{x) obey 

I - 9^ + /i^ + ^ { <Pl{x) + G{x, X) } } Ux) = . (34) 

The ensemble average of all the bilinears of and are determined by the initial ensemble 
average. For our classical evolution we have 

( «fc«fc' )o = ( )o = 27rnfc(/5o) 5{k - k') . (35) 

where, again, rifc(/5o) = l/{(3QUJk), with = a/F + /j?. 

We also need to specify the initial values of fk{x) and fk{x) to complete the calculation. 
These are determined from the initial averages given in Eq. (|TB|) and (00)^, (7i"7r)^, and 
(07r)^ at t = 0. We find: 

/fc(0) = , A(0) = -zujk fk{0) . (36) 

To summarize, for the translationally invariant case when (f)c{x,t) is independent of x, 
the equations that need to be solved simultaneously are: 

[d? + Xiit) ] 0c(t) = , [d^ + k' + X2{t) ] m = , 

Xi{t) =f^' + + y § 2n,(/3o) \Mt)\' , 



— oo 

+ 00 



X2{t) = ^2 + :^^l^t) + '-^ / ^2nM) m)\' . (37) 
2 2 2n 
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The ensemble average of the fields is given by: 

^ -2nM)m)\'. (38) 

We remark here that, with the replacement 2nk + I ^ 2nk, the evolution equations for the 
Hartree approximation in quantum mechanics is the same as the classical ones. In addition, 
in the quantum case, in thermal equilibrium, Uk is given by a Bose-Einstein occupation 
distribution. Quantum mechanics also affects the initial conditions of the mode functions, 
whose Wronskian, in the quantum case, must yield ih. 

For completeness we mention here that in the quantum case, the effective action for the 
Hartree approximation is 

S,s^ J da;£,i(<^,X2) + ^Tr[ln(n + X2)], (39) 
where Ccii<P, X2) is the classical action written in terms of the auxiliary field X2, given by: 

>Ccl(0, X2) = H^) [ □ + X2{X) ] 0(X) + ^ <l>\x) + ^ 

Here the 3 in the Hartree approximation occurring in front of A/2 would become 1 + ;^ when 
there are N fields, showing that the Hartree reduces to the large- A?" approximation at large 
N. 

We notice that unhke the case of the large- A?" expansion, the classical Lagrangian still has 
a quartic term which has to be treated as an external classical field when doing the remaining 
quadratic path integral over to obtain the above determinant term in the effective action. 
Thus, when we take expectation values of the Lagrangian to obtain the energy momentum 
tensor, the quartic term only gets a classical contribution and does not give fluctuation 
contributions. In taking the expectation value in the Hartree approximation one treats x ^ 
well as the quartic term in classically. 

In terms of the mode functions, the energy density is given by: 

e^l [{dt<l>r + {d.<l>r + X2<l>l]-j^ 

+ ^2nM) [\m\' + {k' + X2)\hn (40) 



fJ'^ X2{x) 



xiix) 



fJ'^ X2{x) 
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Because classically rikUJk = To represents the equipartition of energy, the classical energy 
density suffers a linear divergence. In the quantum case, the Bose-Einstein distribution 
cures this, and instead one obtains a quadratic divergence from the zero point quantum 
fluctuations, which must be renormalized by a constant counterterm, corresponding to a 
cosmological constant. 



IV. BARE VERTEX APPROXIMATION 



The Hartree approximation truncates what would be otherwise an infinite set of coupled 
equations for the expectation values of products of the classical field by assuming that all 
fluctuations higher than the second one are exactly zero. In quantum mechanical applica- 
tions, this approximation has many nice features such as being derivable from a variational 
approximation and having a well defined positive definite density function at all times. 
However, it suffers from not containing the hard scatterings usually associated with ther- 
malization. It is thus important to find approximations which include scattering and which 
are free from the problems of secularity and positivity as we discussed in the introduction. 

From the generating functional of all the correlation functions one finds that the Green 
functions for both quantum and classical field theory obey a set of coupled integral equations 
relating the exact two point functions to themselves and higher order one-particle irreducible 
vertex functions. In quantum field theory, the easiest way to derive these is by using a 
generating functional for the matrix Green functions of the CTP formalism. The classical 
theory has similar structure to the Green functions found when we tridiagonalize the matrix 
Green functions and write SD equations for the retarded Green functions and Wightman 
functions. The classical limit {h —>■ 0) of the BVA is obtained by using the classical value 
for the free retarded Green functions and Wightman function and also only keeping those 
graphs which are leading at high temperature and small coupling as discussed by Aarts and 
Smit^^ and by Buchmiiller and Jakovacp3|, |2^ . For the classical theory, one can define free 
propagators which replace the quantum ones by replacing commutators by Poisson brackets, 
as discussed by Parisi|^l|]. In appendix |B], we show how to make precise the connection 
between the exact quantum and classical Schwinger-Dyson equations, by comparing the 
CTP formalism with the MSR formalism. We find that certain bare vertices are missing in 
the classical case. 
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The BVA is an approximation in which the exact two point functions are used in the 
equations for )^ and (0(a;)0(a;') but the exact {(j){x)(f){x')x{x") )^ vertex function is 

replaced by the bare one. This approximation was discussed from several points of view 
in our previous paper |T^, where we showed that the BVA is a particular resummation of 



the expansion. For the classical problem we treat here, the exact SD equations are 
obtained from the MSR formalism |TB|. We sketch the details of this in appendix 0. What 



is important about the BVA is that this approximation is energy conserving, is nonsecular 
and, as Kraichnan has shown [|17|, the effective Hamiltonian has positive spectra. 

In condensed matter physics the BVA is used in the electron-phonon problem for the 
description of the normal state of a Fermi liquid. In this case the first vertex correction 
to the self-energy is small compared to the bare vertex and are of the order of nie/Mion 
where rUe is the electron mass and Mjo„ is the ion mass. This result is known as the Migdal 
theorem |P5||. Close to a superconducting instability the theorem is not valid and a partial 



resummation of the vertex corrections is necessary |^6[|. The physical reason for the validity 
of Migdal's theorem is the adiabatic motion of the electronic degrees of freedom compared to 
the ionic degrees of freedom. In quantum field theory it is not yet clear what the domain of 
validity of the BVA will be. Here we are concerned with its validity in the high temperature 
limit. 

A. The quantum Schwinger-Dyson equations in the BVA 

The exact SD equations for field theory, rewritten in terms of an auxiliary field, 



were first derived by Bender, Cooper and Guralnik|T^ and discussed in detail in the CTP 



formalism in Ref. |T6[. We will also derive the SD equation for both the classical and quantum 
cases in the appendix. The CTP formalism and notation is reviewed in appendix ^ For 
the case where (j){x) = ( (j){x) )^ = 0, one has that 

D{x^ x') = Dq{x, x') — / d^Xi d'^X2 Dq{x, Xi)Il{xi^ X2) D{x2, x') , (41) 

Jc Jc 

G{x,x') = Gq{x,x') — I d^Xi d'^X2Go{x^Xi)Tj{xi,X2)G{x2,x') , (42) 

Jc Jc 



The gap equation for x(t) is: 



xit)=fi' + ^Gix,x)/i. (43) 
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Here Dq {x, x') and Gq {x, x') are given by: 

D,\x,x') = -j6c{x,x'), (44) 
G,\x,x') = [n + xit)]5cix,x'). (45) 

Thus Dq{x,x') = —XSc{x,x'). The exact equations for the polarization 11 and the self 
energy S are 

n(x, x') = 2 y ^^i^^2G{x,Xi)T{xi,x' ,X2) G{x2,x) , (46) 
E(x..')=./dx.dx.G(x.x,)r(x,x'.x.)Z)(x...). (47) 

In the BVA we make the further approximation that the exact one-particle irreducible x00 
vertex function F is replaced by its lowest order value. 

T{xi,x', X2) = 6{xi — x')6{x2 — x') . (48) 

so that we have in the quantum BVA 

U{x,x') = ^G{x,x')G{x',x) , E{x,x') =iG{x,x')D{x',x) . (49) 

To remove the tadpole contributions we write 

D{x, x') = Dq{x, x) + D{x, x') = —A Sc{x, x) + -D(x, x') , (50) 

which we put into Eq. (^Tj), to find an integral equation for D{x,x'): 

D{x, x') = -X^ n(x, + d^x" n(a;, x") D{x", x') . (51) 

In addition, we put 

S(x, x') = A [G{x, x)/i] Sc{x, x') + S(x, x') , (52) 

where 

E(x, x') = i G(x, x') Z)(x', x) . (53) 
Now by multiplying Eq. (^) by Gq ^(x, x') gives a differential- integral equation for G(x, x'): 
[ □ + x{t) ] G{x. x') = Sc{x, x') - [ d^x" S(x, x") G(x", x') , (54) 

14 



c 



Putting (ID into (|5|) gives: 

[ □ + x{t) ] Gix, x') = 5c{x, x') - J d^x" S(x, x") Gix", x') , (55) 
where we define x(t) by 

3A 

m = x{t) + A G{x, x)l% = + Y G{x, x)/i , (56) 
which is the Hartree approximation effective mass. In summary the full set of equations are 

[ □ + x{t) ] G{x, x') = 5cix, x') - d^x" x") G{x" , x') , (57) 

D{x, x') = -A^ n(x, ^') + ^ d^x" n(x, x") D{x", x') . (58) 

f:{x,x') =iG{x,x')D{x',x) . (59) 

U{x,x') = ^G{x,x')G{x',x) , (60) 

x(t)=/i' + yG(a;,x)/z. (61) 

Note the factor of 3/2 in Eq. (|l|). We can now redefine Gq{x,x') using x(t) rather than 
x(t). We call this propagator Go{x,x'). So let 

G,\x,x') = [a + xit)]6cix,x'). (62) 

Then Eq. (|57|) can be written as an integral equation for G{x, x') using Gq{x, x') and S(a;, x'): 

G{x, x') = Go{x, x') — / d^Xi / d'^X2Go{x,Xi)'E{xi,X2) G{x2,x') , (63) 

Jc Jc 

which is completely equivalent to the original equations, except that now we have explicitly 
removed the delta-function term in D{x,x') and included it in the definition of xit)- 

B. Two-particle irreducible eflfective action 

As we have discussed in a previous paper, the BVA can also be obtained by keeping the 
two loop graph in the effective action for the generating functional for the 2-PI graphs |]10 
Namely, the effective action is the twice-Legendre transformed generating functional: 



r[$, g] = 5eiass[$] + ^Tr{ \n[g-']}+ '-TT{G-'m g-i} + r^i^^, g\ . (64) 
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where 

For our problem $ consists of both and x ci^nd is also a matrix in CTP space. Sd is the 
classical Lagrangian written in terms of both (p and x- To obtain the BVA one just keeps the 
two loop graph made from three $ propagators. For the case (0) 7^ 0, the $ propagator will 
be non diagonal in (p-x space. The quantity r2[$,^] has a simple graphical interpretation 
in terms of all the 2-PI vacuum graphs using vertices from the interaction term. When 
{(f))^ = 0, one obtains 

Tlx, g^, V] = 5d,3sM + ^Tr{ In [p-^ ]} 

+ ^Tr{ In [ g-' ]} + ^Tr{D-i V + G"^ g^-2} + [g^, V] . (66) 

where g = {g^p, V} and 

T2[g^,T>] = ~ j dxi j (lx2T>{xi,X2)g<f>{xi,X2)g4,{x2,xi) . (67) 

It is often useful to put this generating functional in the Baym form by defining S = 
g^^ — G^^ to obtain: 

= 5ciass[<f] + '-Ti{\n[g-^]} - '-Ti{T.g} + V2[<^M . (68) 

Eq. ( |68[ ) allows us to relate the graphs in the self energy S to those kept in the effective 
action T2 by means of 

- S = — 

2 5g ■ 

C. Classical limit 

As we discuss in the appendices, the SD equations for quantum and classical evolutions 
are similar in structure, the only difference being that the classical evolution contains fewer 
vertices. The interaction part of the CTP Lagrangian in the original matrix bases is given 
by 

^i = \ - (69) 
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Rotating to the advanced-retarded basis by the transformation 



+ O ^0 ) 



h 

2 

h 

2 



X- ~ X ~ 2 ' 



the Lagrangian becomes 



SO that the bare vertices of the rotated Lagrangian are: 



(70) 



(71) 



7x<^-0 = 7</.x-^ = • • • 4 perms • • • = g , 



In the classical limit the vertices of the last line which have an extra factor of h"^ are missing. 
The self energy in the BVA is given by 



(72) 
(73) 



Sab = lajk G^'- D^"" Jimb 



(74) 



In the above j,l are summed over (f) and vr^ whereas /c,m are summed over x ^ind vr^. 
Similarly the vacuum polarization in the BVA is given by 

1 



Haft = - lajk G^'- G^"" 7zmb 



(75) 



where in the above j, I, k, m are summed over and n^. The inverse Green's functions are 
directly related to the self energies as discussed in Appendix 



V. THE EFFECTIVE ACTION IN THE SINGLE FIELD FORMULATION 

In the single field formulation, we consider an action given by: 

T[cj), G] = ^eiassM + ^Tr{ In [G-i ]} - ^Tr{S G} + T^lcj), G] . (76) 

where now S'dass is the action with the usual 0^ self-interactions. The Hartree approximation 
keeps the two-loop contribution to r2[0, ^] given by: 

3A f 

^2[(P, G]Hartree = ^ / dx [ G{x, x) ]^ . (77) 
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It is easy to go beyond Hartree by using the three-loop contribution: 

r2[0, G] three-loop = J J dXi dX2 [ , X2) (78) 

This approximation was discussed in detail by Calzetta and Hu[^], and recently utilized by 



Berges and Cox|jTj] to discuss equilibration in 1 + 1 dimensional quantum field theory. The 
approximation, which we denote as 2-PI, is lower order than the BVA approximation in that 
it can be obtained from the BVA by keeping only the first term in the equation for D{x, x') 
as a power series in A. That is, if we approximate 

D{x, x') ^ -A^ n(x, x') = -A^ G{x, x') G{x', x) , (79) 

in the BVA, then 

E{x,x') ^ — [G{x,x')f . (80) 

So ( p5|) reduces to: 

\U^x{t)\G{x,x') = bc{x,x')-— [ d^x"[G{x,x")]^G{x",x') . (81) 

2 Jc 

We can rewrite this as: 



|n + + ^Gy{x,x)/i I G>(a;, a;') 



y I dx"{ I dt"G>{x',x")Gl{x,x")+ I dt"G>ix",x')Gl{x,x") 



dt" G> {x", x') Gl (x", x) I , (82) 



which agrees with Eqs. (9) and (10) of Berges and CoxljT^ apart from a symmetry factor 



of 1/3 on the right hand side of Eq. (|82D . For the classical case, we need to only include in 
n>(x,x') and Il^{x,x') the subset of vertices that come from the classical MSR formalism 
(see appendix 0). What we find is that the results for the 2-PI approximation are quite 
similar to those found for the BVA. 



VI. RESULTS 

The present codes for evolving the SD equations use spectral numerical methods based 
on Chebyshev polynomial expansions, as explained inp8|, p9[ . Thus our calculations are not 
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lattice calculations in coordinate space using periodic boundary conditions, but are carried 
out entirely in momentum space. They take advantage of the global spectral convergence 
character of the Chebyshev polynomial expansion, which allow for an exact continuum limit. 
The calculation of the different modes decouple and an implementation via the message 
passing interface [MPI] on a parallel computer is straightforward. 

We choose our initial parameters to match those given by Aarts, et We have taken 

A = 1/3 and /i = 1, which are the values which make the classical problem parameter free. 
We set A = vr/a = Att^ = An. This makes a = 1/4. The physical size of the system is then 
fixed by the value of from the relation, L = Na = ttN/A = N/A. 

Also, to make contact with the simulations of Ref. [Q, we choose our initial density 
matrix to be described by Eq. ( pJ]) with the mass parameter set equal to the bare mass 
The value of jSo is taken to agree with the lattice coordinate space calculation of 



Ref. 0, l/Po = Tq = 5.03891094, so that for the coordinate space Monte Carlo calculations, 
( 0) )p = 2.5 exactly. (This is for convenience only.) 
For the Hartree and BVA approximations, we take: 



/fc(0) = 1/V2^k , fk{0) = -luJk /fc(0) = ~i , (83) 

where ujk = a/A;^ + xlO); with 

x(O)=/.^ + y(0^(O))^, (84) 

and /(A) is given by Eq. (pO|). 

In Fig. |1|, we show (l)\{t) = { (^^{t) ) clS cl function of time, using the coordinate space 
lattice calculations of Eq. ( P^D together with the various approximations presented in this 
paper. The BVA, Hartree, large-A^, and coordinate space Monte Carlo codes use quite 
modest values of A^ = 128 or 256, as the calculation is very insensitive to this parameter. 
The momentum space Monte Carlo code uses somewhat larger values of A^ to obtain high 
accuracy. 

We see that that the initial value different for the lattice coordinate and momen- 



tum space approaches, as explained in Section |II D| . For fixed values of the parameters (3q 
and A, the value of 0^1 (0) different in the two formalisms. One can of course force ( 0^(0) )^ 
to be the same, but then the mismatch for finite cutoff would appear somewhere else. Or 
one may try to correct this problem by shifting the lattice result such that it matches the 
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continuum approximations at t = 0. We show the results of such a shift in Fig. |^. However, 
such a naive solution has no fundamental support and may be misleading. 

In order to compare directly our approximation methods to the exact Monte Carlo results, 
we have solved the classical Monte Carlo numerically in momentum space, as explained in 
Section [11 D| . In Fig. ^ we illustrate these lattice momentum space results for and 
compare them with the approximations presented here. The BVA result fares remarkably 
well and cures the oscillation problems of the Hartree approximation at large times. It does 
not display the re-amplifications (secularity) problems manifested by the result of Aarts et 
al, where the next-to-leading order evolution is obtained by truncating the effective action up 
to four-point couplings]^. Similar behavior is noted for the time-evolution of depicted 
in Fig. |, and the effective temperature Tgfj = (7r^(t) ), shown in Fig. |^. 

In Fig. ^ we compare the BVA result with the result of a calculation for which we 
truncate the D{x,x') equation at leading order, which makes it similar to the procedure 



advanced by Berges and Cox|14]. The results are very close indeed. The reason seems to 
be that the Berges approximation is expected to give reasonable results at short times and 
the differences should build up as the time-evolution progresses. However, at later times 
the thermalization mechanism takes over and virtually obliterates the differences. Since the 
classical problem scales, we have only computed results for /i = 1 and A = 1/3. In the 
quantum case, a second dimensionfull parameter (h) enters, so that the problem no longer 
scales, and one will have to treat the case of large coupling constants (A = 7.3 for the linear 
sigma model) and high temperatures. Then the intermediate domain, between the regime 
when the Berges' leading order approximation for D{x,x') holds, and the regime when the 
thermalization blurs the differences, will become more important. 



VII. CONCLUSIONS 



In this paper we have compared exact numerical simulations of classical scalar 1-1-1 di- 
mensional field theory with several Schwinger-Dyson truncation schemes. We have focused 
on approximations that do not suffer from obvious problems, such as secularity and negative 
probability, by obtaining them from self-consistent approximations to the effective action 
which have been shown by Kraichnan to correspond to physically realizable systems. In 
particular, we have shown that a resummation of the 1/A^ expansion that we studied earlier 
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FIG. 1: The lattice coordinate space Monte Carlo calculation of 'p'^iit) = {<p'^{t) )g, compared with 
the leading order large- A^, Hartree, and BVA approximations. 




FIG. 2: The shifted lattice Monte Carlo calculation of <p'^i{t) = { (p'^it) )q, compared with the leading 
order large-A^, Hartree, and BVA approximations. Here, a constant factor is subtracted from the 
lattice coordinate space Monte Carlo results of Fig. |l[ 
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FIG. 3: The lattice momentum space Monte Carlo calculation of ^ i^'^i^) )o' compared with 

the leading order large- A^, Hartree, and BVA approximations. 
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FIG. 4: The lattice momentum space Monte Carlo calculation of compared with the leading 
order large- A^, Hartree, and BVA approximations. 
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FIG. 5: The effective temperature, Tgff, as defined in tlie text, for tlie lattice momentum space 
Monte Carlo calculation and the BVA result. 




FIG. 6: (p'^iiO) = )q for the BVA approximation and the approximation of Berges and Cox, 

for = 1, A = 1/3. Note the close agreement. 
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in a quantum mechanics model |T^, the bare vertex approximation [BVA], cures the late 
time oscillation problem of the Hartree approximation, and leads to results that quantita- 
tively agree with exact numerical simulations, even at late times. Since the classical field 
theory, averaged over an initial Maxwell-Boltzmann distribution, is the high temperature 
limit of the quantum field theory, averaged over an initial Bose-Einstein distribution, our 
results imply that the BVA when applied to the quantum evolution meets the minimum 
requirement of being valid in the high temperature domain. 

Our results give us confidence that this approximation will be useful in future studies of 
quantum phase transitions in the 0(4) model. The fact that in homogeneous situations this 
approximation leads to thermalization, will allow us to study the rate of equilibration versus 
expansion rate for an expanding plasma undergoing a phase transition. We will then be 
able to see whether some of the interesting phenomena occurring during the early stages of 
a chiral phase transition will survive the hard scatterings present in the BVA approximation. 
At the same time this study was done, a complementary study of this approximation in the 
context of the quantum field theory in 1+1 dimension was performed with results not that 



different from the classical case|3^. Again, in the quantum simulations, the oscillations 



present in the Hartree approximation were found to damp, and equilibration occurred. 
APPENDIX A: CTP FORMALISM 

For the quantum field theory problem, the SD equations for initial value problems can 
be obtained directly from the CTP path integral for the generating functional. The 
functional differential equation for the generating functional is exactly the same as that of 
ordinary quantum field theory, with the exception that the ©-functions used in defining 
time-ordered products are defined on the closed-time path contour, rather than on the real 
line. This leads to a 2 x 2 matrix formulation for the Green functions ||19||. It is useful to 



have three different forms of the CTP formalism. The first is the original matrix formulation 
of Schwinger, the second uses closed contour 9-functions, and the third is a rewrite of the 
matrix formalism in terms of Wightman functions and advanced and retarded propagators. 
The generating functional for initial value problem Green's functions is 

Z[J+, J-] = e*^I''"l = J lldipaexpt{SM + J^ip''}{iput\p\ip2,t), (Al) 
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where ( (/?!, i| p |</?2, is the density matrix defining the initial state. We use 
notation 



for a = 1, 2. 



On this matrix space there is an indefinite metric 



Cab = diag (+1,-1) = c' 



SO that, for example. 



J" Cab = J+ - J- 



Prom the path integral we get the following matrix Green's function: 



SJa{t) 6Jh{t') 

We can also write this using CTP 0-functions as 



j=o 



G{x,x') = Tr{ Tc[4>{x)(l>{x')]} = Qc{t,t') G>{x,x') + Qlit' ,t) G^{x,x'] 

x') G {x,x') 

where the CTP ©-functions are defined by: 



1 iort> t' on C, 
for t < t' on C, 
on the time contour C. In matrix form, we have: 



which leads to 



^e(t - 1') 
^ 1 e{t' - 1)^ 



e{t' - 1) 1 

Q{t - 1') 



G++{x, x') = Q{t - t') G>{x, x') + e(t' - t) G<{x, x') , 
G_+{x,x') ^G>{x,x') 
G+_{x,x') = G^{x,x') 

G^4x, x') = e{t' - t) G>{x, x') + e{t - t') G<{x, x') , 
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where 



G>{x,x')/i = Tr[p0(x)0(x')] , G<{x,x')/i = Tr[p0(x')0(x) ] , 



(A9) 



with <j){x) and p quantum operators. Here, G>{x^x') satisfy the symmetry relations 



Gy{x,x')/i = Gy{x',x)/i = G<{x,x')/i 



(AlO) 



If we use the Schwinger-Keldysh basis then we define the 000 vertex as: 
|/i„bc0"0'0% Kbc^±l, ifa = 6 = c = ±l. 
The two point function SD equation can be written as 



(All) 



Gab /^ab r^ac v r^dh 



(A12) 



To make a direct relation with classical correlation functions, it is convenient to go to another 
basis, given by 



01 

,02, 



(0+ + 0_)/2 

0+ - 0- 



R 



(A13) 



with 



R 



1/2 1/2 



(A14) 



To make contact with the MSR formahsm discussed below, we will choose 0i = and 
02 = /iTT,^, so that 



0± = ± \^'l>- 



(A15) 



The usual matrix for the propagator 



G(x - x') 



Gj^j^{x — x') Gj^Ax — x' 



G_^[x-x') G__{x-x') 



(A16) 



gets transformed into 



G(x -x')^R G(x - x') IV = 



iF[x — x') Gr(x — x' 
Gdx-x') 



(A17) 
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with 

F{x -x') = ^{ 0(x) + ) , 

Gr(x-x') =Ga{x' -x) =i{[(f){x),(j){x')])e{t-t') . (A18) 

The last equation can be used to define the spectral function: 

a{x -x') = i{ [(j){x), (f){x') ] ) = Gr{x - x') - Ga{x - x') . (A19) 
In the advanced-retarded basis we then get the integral equations: 

F{x, x') = Fo{x, ^')- Jj dx"dx"' { Goa{x, x") Sf(x", x'") Gr{x"', x') 

+ Fo{x, x") Sr(x", x'") Gr(x"', x') + Goa{x, x") SA(a;", x'") F{x"', x')} 
Ga{x, x') = GoAix, x') — j j dx"dx'" Gqa{x, x") Tja{x'\ x'") Ga{x"', x) 

GrIx, x') = GoR(a;, x') - ! ! dx'dx'" Gor(x, x") Sr(x", x'") Gr(x'", x'). (A20) 



The SD equations for our constrained x fi^ld problem can be obtained by considering first 
two propagating fields, and then taking the composite field limit for the second field, so 
the bare propagator for the second field is replaced by a delta function (see Cooper and 
Haymaker Since the second approach is more transparent in the MSR formalism we 



will follow that here. Thus we start with the Lagrangian 

^ = \{ + id,x? -m'^'-M\^-gx^'}-J(p- (A21) 

This leads to the equations of motion in the presence of external sources: 

[n+m^]^{x) = -g x{x) 0(x) - j (x) . (A22) 

[ □ + ] xix) = -g 02(x)/2 - S{x) . (A23) 

After we determine the SD equations for the two field problem, we recover the original 
equations for the single 0^ field theory by taking the composite limit where: 

xix) = + 

Don{x~y)^-X5{x-y). (A24) 
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We want to convert the Schwinger-Keldysh Lagrangian to the advanced retarded one. Start- 
ing from: 

L = \ {{d.ct^^f - id,<P^f + id,x+f - {d,X-f - (0^ - - {xl - x'-) 

- 9 ix+ 0+ - X- 0-) } - - 0-) - S{x+ - X-) ■ (A25) 

we obtain after changing variables to 0, tt^ 



-LcTP = vr4 □ + ] + vr^ [ □ + ] X + TT^ X + I vr^ 0' 



+ ^9^<t>^x- H(t>- JxX- Jir^ T^<p - J^x ■ (A26) 

We will find that we get exactly the same Lagrangian in the MSR formalism except that 
the term proportional to ff' is missing. 

APPENDIX B: MSR FORMALISM 

In the paper of Martin-Siggia-Rose [MSR] |jl8| , an operator formalism was developed which 
allowed them to find a generating functional for both the retarded and Wightman functions 
for first order classical field equations of the type: 

x{r,t) = A[x{r,t)] (Bl) 

where y4[a;(r, t)] is a local polynomial in the classical field x{r,t). In the work of MSR, 
A[x(r, t)] could contain dissipative terms as will as prescribed noise terms. The formalism 
presented in MSR is first order in time derivatives and not apparently covariant. We have 



recently shown pO| that there is a covariant subset of the MSR equations in terms of which 
all the MSR Green's functions can be recovered. This subset can be derived from a second 
order Lagrangian formulation which can be related to the h —>■ limit of the CTP formalism 
of Schwinger and Keldysh as mentioned above. 

1. First order formalism 

First let us review the first order formalism. For the statistical classical field evolutions 
of interacting classical fields (pa, a = 1,2, . . . , N then x{r,t) is the 2N component field 
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consisting of (pa and the canonical momentum TTa — (j)a- 



X 



(B2) 



7r„, 



If, for example, we restrict ourselves to cubic interactions, then the vector A is of the form: 

Ai = Ci{r, t) + dijXj{r, t) + eijkXj{r, t)xk{r, t)/2 (B3) 

where i — 1, . . . ,4. In the MSR formahsm one then introduces the operator 

x(r, t) = —-r-, r 

^ ' ^ 5x{r,t) 

such that the commutation rule [x{r,t),x{r' ,t)] = 6{r — r') is true. Then if we define an 
operator Hamiltonian via: 

H{t)^ j dr' Xi{r' ,t) Ai{r' ,t) , (B4) 
then the equations of motion can be written in the compact form 



x{T,t) = [x{T,t),H{t)] 

For the commutator to be true at all times, one needs that x satisfies 

dx(r, t) 



[x{T,t),H{t)] 



(B5) 



(B6) 



for consistency. Therefore x{r,t) is a functional of x(r, 0) and x(r, 0). The formal solution 
to these equations is given by 



x{t) = C/-^(i,0)x(0)C/(i,0) 
x{t) = C/-^(i,0)x(0) U{t,Q) 



(B7) 



where 



C/-i(t,to) = T|exp 
C/(i,io) = T*|exp 



H{t') dt' 



to 



+ f H{t')dt' 

•'to 



Here T corresponds to the usual time ordered product operation. The meaning of the 
expectation value {x{t) x{t')) is as follows. Given an initial probability function P[x(0)], 
then 



( x{t) x{t') ) = / da;(0) x[t, x{0)] x[t', x{0),x{0)] P[x{0)] . 



(B8) 
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Wherever i;(0) appears, it is replaced by 

x(0) 



dx(0) 

and it acts on everything to the right. This definition of the extended averaging procedure 
has three important properties. 

1. The average of a product of x's agrees with the conventional definition. 

2. The time dependence of {x(t) x{t') ) is consistent with the above equations of motion. 

3. The expectation value of a product of x and x which has an x to the left vanishes if 
P[a;(0)] goes to zero fast enough at large 

The last property is crucial for the tridiagonal form of the Green functions and follows from 
the fact that 

/ [x"P(x) ] da; = , if: lim x"P[a;] = . 

J-oo dx \x\^oo 

Thus in particular {x{t')x{t) ) = 0. 

The meaning of the hatted operators is understood in terms of the response to the system 
to an external source, as is discussed in Refs.|]l8[ and |^ . 



2. Second order formalism 

Let us first write the equations for the two coupled fields x and in first order form. We 
have 

§ = ^=[V^ + M^]x + 9<PV2. (B9) 

The operator Hamiltonian which generates these equations is 

^^MSR = j dx 1^ TT,^ + x% 

-n^[[\/' + m^]<P + gX^}-n^{['^' + M']x + g<pV2]\ (BIO) 
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Following the arguments of Ref . pO| , and using the commutation relations, we find that the 



independent covariant second order equations (adding sources) are 

[ □ + ] + 5( X </> = J</, , 

[□ + M2]7r^ + (707r^ = j^^. (Bll) 

Here 7r(^(x, 0) is treated as the operator 6/6(j){x, 0) when one averages over the initial prob- 
ability function in phase space. These equations are derivable from the Lagrangian density: 

which is identical to ( |A26| ), except that the terms proportional to fi^ are missing. This is 
analogous to what we showed for the quantum mechanics example in Ref. 0]. The extra 
bare vertices of the quantum field theory calculation are obtained from the term: 

%9^l^x- (B13) 



3. The classical SD equations 

The way one derives the SD equations from the action is identical for both classical 
and quantum field theory. Thus we obtain the same structure classically, but there are 
fewer vertex functions. We now derive the SD equations for a generic cubic self-interacting 
field theory whether classical or quantum. For interacting scalar fields, we introduce the 
column vector $q, composed of $J = (j) and $^ = txi where i = 1,2,... ,N. We also need 
the metric ga/3 which is just cr^^ as far as connecting the $^ and sectors. Then we can 
write generically for cubic interactions: 

£ = i D^l^ $^ + ^ 7,^, $° $^ <I>P - J, (B14) 

where 

DoU^ -y) = gap[a + m'] 6{x - y) . (B15) 
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The generating functional is formally 



W[J] 



Z[J] = (T{exp[J,$"]}) = e 



Defining the "classical" field $° and the connected 2-point function ^2(1, 2) via 



$"(1) 



5W 



one has that in the presence of sources 



Woiix 



SM2) ' 



6U2) 



+ $^(1) W'^'il, 2) + $"(1) Wf^Pil, 2) 



Jail) , 

(B16) 



To obtain the SD equation for the inverse two point function we first use the connection 



between the connected 3-point function and the 1-PI vertex function||10|| obtained by first 
Legendre transforming from J to $ and using the chain rule 



SJJl] 



d2iy(l,2 



,«/3 



/3 



(B17) 



We find: 



(r{$.(i)<i.^(2) $,(3)}) 



5W^pil,2) 
6JP{3) 

d4 d5 d6 W^^,{1, 4) Wf,pi2, 5) H^,,,(3, 6) r-'^''''(4, 5, 6) 



where the 1-PI three point function is defined by 

5W-Ul,2] 



ra/3p(l, 2, 3) 



0/3 V-^) ■ 



5$^(3) 



(B18) 



(B19) 



The usual SD equation is obtained by multiplying Eq. ( |B16| ) on the right by ^ to obtain 



WJ{1, 2) = D-Ul, 2) + S,^(l, 2) + 2 7,^, $"(1) 6(1 - 2) , 



(B20) 



with 



S„^(l,2) = j d3d47„,,W^''\l,3)l^'^'^(l,4)r(3,4,2; 



(B21) 
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and where 



r,.,(l, 2, 3) = 2 7„.p 5(1 -3)5(1-2) + -j^^ ■ (B22) 

However S is a proper self energy and can be related^] to the effective action T2[G] for the 
2-PI graphs of the theory via 

SO that it is just a function of the bare vertices and the full Green functions. Using 



5$p(3) J 5^-^(5,6) 6^p{3) 

we find that 

r„^,(l,2,3) = 27„^,5(l-2)5(l-3) 

- y"d4d5d6d7r„,,(l,4,5)G^^(4,6)G^\5,7)i/^A;/3.(6,7;2,3). (B25) 

The scattering kernel H is defined by 

^^^^-if„,.,(l,2,5,6). 

Self consistent approximations are determined by keeping a certain class of graphs in 
r2[G'], the sum of all 2-PI graphs made from bare vertices and full propagators. For cubic 
interactions, the BVA is obtained by keeping the graph 



j dl d2 7^,,G^'(1, 2) G^™(1, 2) G^"(l, 2) 7^, 



which then leads to the self energy being the one loop diagram, and the scattering kernel 
being single particle exchange. By using the variational definitions of S and H one is 
guaranteed an internally consistent approximation. In the BVA, the integral equation for 
the correlation function is 

^,^(1,2) =G'oafc(l,2)+ j d3MGoa^{l,3)l^,kG''{3,A)G''^{3,4)lln.nG^b{^,2). (B26) 

Expanding (P26|) , we can write this in terms of classical and quantum contributions, as 
shown in ||2^. We identify the self energy in this approximation as 

S,„(3, 4) = 7,j, ^^'(3, 4) G^"^(3, 4) 7;^„ . (B27) 



33 



ACKNOWLEDGMENTS 

Present calculations are made possible by grants of time on the parallel computers of 
the Mathematics and Computer Science Division, Argonne National Laboratory. We thank 
Harvey Rose for explaining the MSR formalism and for providing us a copy of his unpublished 
thesis. JFD and BM would like to thank LANL for hospitality during part of this work. The 
work of BM was supported by the U.S. Department of Energy, Nuclear Physics Division, 
under contract No. W-31-109-ENG-38. 



[1] A. K. Kerman and S. E. Koonin, Ann. Phys. 100 (1976) 332; R. Jackiw and A. K. Kerman, 
Phys. Lett. A 71 (1979) 158; A. H. Guth and S.-Y. Pi, Phys. Rev. D 32 (1985) 1899; F. Cooper, 
S.-Y. Pi and P. Standoff, Phys. Rev. D 34 (1986) 3831; S.-Y. Pi and M. Samiullah, Phys. 
Rev. D 36 (1987) 3128. 

[2] D. Boyanovsky and H. J. de Vega, Phys. Rev. D 47 (1993) 2343; D. Boyanovsky, H. J. de Vega, 
R. Holman, D.-S. Lee, A. Singh, Phys. Rev. D 51(1995) 4419; D. Boyanovsky H. J. de Vega, 
R. Hohnan, J. Salgado, Phys. Rev. D 54 (1996) 7570; D. Boyanovsky, D. Cormier, 
H. J. de Vega, R. Holman, A. Singh, M. Srednicki, Phys. Rev. D 56 (1997) 1939; D. Boy- 
anovsky, M. DAttanasio, H. J. de Vega, R. Holman and D. S. Lee, Phys. Rev. D 52 (1995) 
6805; D. Vautherin and T. Matsui, Phys. Rev. D 55 (1997) 4492; D. Boyanovsky, H. J. de Vega, 
R. Holman and J. Salgado, Phys. Rev. D 57 (1998) 7388. 

[3] F. Cooper and E. Mottola, Phys. Rev. D 36 (1987) 3114; F. Cooper, Y. Kluger, E. Mottola 
and J. P. Paz, Phys. Rev. D 51 (1995) 2377; Y. Kluger, F. Cooper, E. Mottola, J. P. Paz, 
A. Kovner, Nucl. Phys. A590 (1995) 581c; M. A. Lampert, J. F. Dawson and F. Cooper, Phys. 
Rev. D 54 (1996) 2213; F. Cooper, Y. Kluger, and E. Mottola, Phys. Rev. C 54 (1996) 3298. 

[4] A. Chodos, F. Cooper, W-j. Mao, H. Minakata, A. Singh, Phys. Rev. D61 (2000) 045011 
(|hep-ph/9909296|| ; A. Chodos, F. Cooper, W-j. Mao and A. Singh, Phys. Rev. D63 (2001) 
096010. 



[5] C. Wetterich, Phys. Rev. Lett. 78 (1997) 3598 [|liep-th/9612206|] ; L. Bettencourt and C. Wet- 
terich, Phys. Lett. B430 (1998) 140 [|hep-ph/9712429| ; L. Bettencourt and C. Wetterich, 
||hep-ph/9805360t ; G. F. Bonini and C. Wetterich, Phys. Rev. D 60 (1999) 105026 



33 



ph/9907533]; 

G. Aarts, G. F. Bonini and C. Wetterich, Phys. Rev. D 63 025012 (2001) [hep-ph/007357]. 
F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. Paz, and P. Anderson, Phys. Rev. D 50 
(1994) 2848 f|hep-ph/9405352|| ; F. Cooper, J. F. Dawson, S. Habib, Y. Kluger, D. Meredith 
and H. Shepard, Physica D 83 (1995) 74. 

B. Mihaila, T. Athan, F. Cooper, J. F. Dawson, and S. Habib, Phys. Rev. D 62 (2000) 125015 
[|hep-ph/0003T05 |. 



J. M. Luttinger and J. C. Ward, Phys. Rev. 118 (1960) 1417; G. Baym, Phys. Rev. 127 
(1962) 1391; B. Vanderheyden and G. Baym. To appear in Progress in Nonequilibrium Green's 
functions, M. Bonitz (Ed.), World Scientific, Singapore 2000 [|liep-ph/000"229l|. 



M. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D 10 (1974) 2428. 

C. De Dominicis, J. Math. Phys. 3 (1962) 983; C. De Dominicis and P. C. Martin, J. Math. 

Phys. 5 (1964) 14, ibid, 5 (1964) 31. 

H. D. Dahmen and G. Jona-Lasinio, Nuovo Cim., 52A (1967) 807; ibid, 62A (1969) 889. 
A. N. Vasil'ev and A. K. Kazanskii, Teor. Mat. Fiz., 12 (1972) 352; ibid, 14 (1973) 289. 
J. Berges and J. Cox, "Thermahzation of Quantum Fields from Time- Reversal Invariant Evo- 
lution Equations," |hep-ph/0006T6C| ]. 



C. M. Bender, F. Cooper, and G. S. Guralnik, Annals Phys. 109 (1977) 165. 
B. Mihaila, F. Cooper, and J. F. Dawson Phys. Rev. D 63 (2001) 096003 f|hep-ph/ 000'6254|] . 
R. H. Kraichnan, J. Fluid. Mech. 5 (1959) 497; J. Math. Phys. 2 (1961) 124; J. Math. Phys. 
3 (1962) 205. 

P. C. Martin, E. D. Siggia and H. A. Rose, Phys. Rev. A 8 (1973) 423; H. A. Rose, Harvard 
Thesis (1974) (unpublished). 

J. Schwinger, J. Math. Phys. 2 (1961) 407; P. M. Bakshi and K. T. Mahanthappa, J. Math. 
Phys. 4 (1963) 1; 4 (1963) 12; L. V. Keldysh, Zh. Eksp. Teo. Fiz. 47 (1964) 1515; [Sov. Phys. 
JETP 20 (1965) 1018]; G. Zhou, Z. Su, B. Hao and L. Yu, Phys. Rep. 118 (1985) 1; 
F. Cooper, A. Khare and H. Rose, "Classical limit of time-dependent quantum field theory 
— a Schwinger-Dyson Approach," Phys. Lett. B515 (2001) 463 [|hep-ph/0106Tl^ . 



G. Parisi, Statistical Field Theory, (Addison- Wesley, New York, 1988). 
G. Aarts and J. Smit, Phys. Lett. B393 (1997) 395 [lhep-ph/9610415| . 



W. Buchmiiller and A. Jakovac, Phys. Lett. B407 (1997) 39 ||hep-ph/9705452|] . 



35 



[24 
[25; 
[26 
[27 

[28 
[29 
[30 



W. Buchmiiller and A. Jakovac, Nucl. Phys. B521 (1998) 219 [|hep-th/9712093 
A. B. Migdal, Sov. Phys. JETP 7 (1958) 996. 

J .R. Schrieffer, Theory of Superconductivity, (Benjamin, New York, 1964). 

E. A. Calzetta, and B. L. Hu, Phys. Rev. D 37 (1988) 2878; E. A. Calzetta, B. L. Hu, S. A. 

Ramsey, Phys. Rev. D 61 (2000) 125013. 



B. Mihaila, J. F. Dawson, and F. Cooper, Phys. Rev. D 56 (1997) 5400 |hep-ph/9705354 |. 
B. Mihaila and I. Mihaila, [ |physics/990"T005 |. 

J. Berges, "Controlled nonperturative dynamics of quantum fields out of equilibrium," [ |hep-| 



ph/0105311|] . 



[31] F. Cooper, G. Guralnik, R. W. Haymaker, K. Tamvakis, Phys. Rev. D 20 (1979) 3336. 

[32] W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery, Numerical Receipes in Fortran, 
2nd edition, (Cambridge University Press, 1992), p. 833. 

[33] This choice is not optimal, since the effective mass at t = is the Hartree value. In classical 
1+1 dimensional field theory there is only finite mass renormalization, and using ^u^ in our 
initial density matrix is allowable. However, if we were to study the quantum version of this 
problem, we would have to use the renormalized mass and not the bare mass in the initial 
density matrix. 



36 



